Casimir Torques between Anisotropic Boundaries in Nematic Liquid Crystals 
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Fluctuation-induced interactions between anisotropic objects immersed in a nematic liquid crystal 
are shown to depend on the relative orientation of these objects. The resulting long-range "Casimir" 
torques are explicitely calculated for a simple geometry where elastic effects are absent. Our study 
generalizes previous discussions restricted to the case of isotropic walls, and leads to new proposals 
for experimental tests of Casimir forces and torques in nematics. 
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sions in critical fluids or superfluid s [s] . 




In the last decade, much theoretical attention has 
been paid to "Casimir" forces in structured complex flu- 
ids |-|. The pioneering work of Casimir showed that 
two uncharged conducting plates attract each other in 
vacuum, due to the modification of the electromagnetic 
fluctuations imposed by the plates In complex flu- 
ids, similar interactions should exist between embedded 
objects, as the thermal fluctuations of the medium's elas- 
tic distortions are restrained by the boundary conditions 
imposed by the objects. These interactions are believed 
to act between bounding surfaces or immersed inclu- 

, in liquid crys- 
and also be- 
Nematic liquid 



tals [|7[-|10[, in bilayer membranes 
tween rodlike polyelectrolytes [|lT 
crystals are anisotropic fluids with a quadrupolar long- 
range order. They are considered as good candidates for 
the direct observation of "Casimir" interactions in com- 
plex fluids. However, clear experimental evidences have 
to date been scarce |g] . This is probably due to the weak- 
ness of fluctuation-induced effects when compared to that 
of permanent elastic distortions, which are often present. 

Although nematic liquid crystals are well-known to dis- 
play orientational effects |]18[, no "Casimir" interaction 
directly connected to these orientational properties of ne- 
matics have so far been discussed. Here, we demonstrate 
that thermal fluctuations in nematic liquid crystals can 
induce torques between bounding surfaces The ex- 
istence of "Casimir" torques between objects embedded 
in complex fluids is usually caused by the anisotropy in 
the shape of the objects [ p^|p^Jl7[ |; here we report on a 
more subtle effect occurring between infinite plates with 
translational symmetry. To emphasize the "Casimir" ef- 
fect, we focus on a situation in which no average elastic 
distortion is present: two parallel plates with a surface 
energy favoring an orientation of the local average molec- 
ular alignment (director) perpendicular to the surfaces. 
The ground state is therefore the distorsionless state in 
which the director is everywhere perpendicular to the 
boundaries. A "Casimir" torque can arise due to the 
anisotropy in the rigidity of the surface energy: we as- 
sume that deviations from the preferred normal surface 
orientation is easier in one direction than in the orthogo- 
nal one. (This property can be experimentally obtained 
from a grating surface treated for homeotropic anchor- 



ing [gO[, or by depositing on top of a substrate which is 
conventionally treated to give planar anchoring a very 
thin layer of a material promoting homeotropic anchor- 
ing ||2l|.) In this geometry, we calculate the "Casimir" 
interaction, and show that it depends not only on the dis- 
tance between the two plates, but also on their relative 
orientation. At the end of the paper, we argue that this 
leads to effects easier to measure than the direct force 
between the plates. 

Nematic liquid crystals are liquids of rodlike molecules 
displaying a long-range orientational order. The local 
average molecular axis is represented by a unit vector 
±n, called the "director" . The bulk ground state corre- 
sponds to a uniform director field and the Frank elas- 
ticity describes the free energy associated with gradi- 
ents of the director ||2^. Bounding plates often favor 
some orientation of the director: this phenomenon is 
known as anchoring JTst . The simplest situations cor- 
respond to a preferred orientation normal to the plates 
(homeotropic anchoring) or parallel to the plates (pla- 
nar anchoring). Here we employ a path integral method 
to study the "Casimir" energy of a nematic liquid crys- 
tal confined between two parallel plates at a separation 
d, at which we assume an homeotropic anchoring with 
anisotropic strength as described above. We calculate the 
fluctuation-induced interaction between the two plates 
when the axes of weakest anchoring strength are placed 
at a relative angle 9 (see Fig. 1), and find 
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per unit area of the plates. The angle dependence in this 
interaction shows that the plates experience a long-range 
fluctuation-induced torque, that decays algebraically as 
l/d^, and tends to align the directions of weakest rigidity. 

We start with the broken symmetry configuration of 
lowest energy, in which the director field is aligned along 
the z-axis, and restrict ourselves to small fluctuations 
(5n = n — z of the director around this ground state (see 
Fig. 1). The bulk cost of such a fluctuation can be de- 
scribed by an effective Hamiltonian 
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which corresponds to the usual one-constant (K) approx- 
imation of the Frank elasticity [^,|8| . The anisotropic 
homeotropic anchoring surface energy is accounted for by 
the surface Hamiltonian 
1 



Hs = 



(fr^ ( 5n- ■ W" • Sn- + Sn+ ■ W+ ■ Sn+) 



(3) 

where (5n^(r^) — Sn{r±,0) and (5n+(rx) — Sn{r±,d), 
and W~ and W'*' are positive definite constant tensors 
that entail the anisotropy of the surfaces and their rel- 
ative orientation. We assume that the two surfaces are 
identical in nature and we define W~ = Wmax x"® x~ -I- 

W^miny"(8>y- and W+ = WmaxX+(^x+ w^mi„y+®y+, 

where x* and y^ denote the hard and weak axis on 
plate ±, respectively (see Fig. 1). The eigenvalues Wmax 
and Wniin naturally define two extrapolation lengths 
= K/Wynax and Amax = ^/Wnin |l8|l ■ Wc assumc 
extreme anisotropy, namely Amin ■C Amax- 




FIG. 1. Sketch of a fluctuation of the bulk director fleld 
n(r) (thin line) around its mean value z (bold line). The 
two anisotropic plates impose a homeotropic anchoring (i.e., 
the preferred surface value of n is z) characterized by an 
anisotropic anchoring elasticity: deviations from z are easier 
in the "soft" directions than in the "hard" directions x^. 
In this geometry, the fluctuation-induced "Casimir" interac- 
tion between the plates depends both on their separation d 
and their relative angle 9. 

"Casimir" effects arise from the quantization of the 
fluctuation modes by the boundaries (essentially from the 



suppression of modes of wavevector smaller than 27r/d). 
The actual effect of the boundaries on a fluctuation mode 
of wavevector q ~ 2Tr/d is a function of the product qX, 
where A is the extrapolation length corresponding to the 
considered polarization of the fluctuation. Depending on 
the relative values of Amin, -^max and d, three different 
regimes occur: (i) For Amin *C Amax *C d the constraints 
are effectively hard in both directions and the anisotropy 
is washed out at leading order, (ii) For Amin ^ d ^ Amax 
the director field is subject to a hard constraint in the 
X direction while it is virtually free in the (perpendic- 
ular) y direction, (iii) Finally, for d ^ Amin ^ Amax 
both directions allows almost free fluctuations and the 
anisotropy is lost again at leading order. To emphasize 
the effect of the anisotropy, we thus focus on the case 

where Amin <C (i < Amax- 

Following Ref. 0, we now proceed to calculate the par- 
tition function for the director fluctuations, using TC + TLs 
as the total Hamiltonian. The quantum mechanical de- 
scription of the partition function Q , which treats z as 
an imaginary time variable, helps us gain a useful insight 
into the meaning of the boundary conditions. From the 
"imaginary time action" 7i, one can define the "momen- 
tum" conjugate to a component of Sn as Psn = K dzSn. 
Then, one can go from coordinate space to momentum 
space. In particular, at the boundaries one finds out that 
Psn is quadratically confined by the tensor W^^ and thus 
acts opposite as compared to 6n. In other words, there 
is an uncertainty principle relating the fluctuations of 6n 
and d^Sn as A{Sn)A{dJn) ~ ksT/K. In hght of this, 
one can argue that the boundary condition of nearly free 
fluctuations in the soft direction is asymptotically equiv- 
alent to setting dzSn — 0. Note that this is equivalent to 
assuming that the directors cannot bear torques due to 
the freedom of rotation in the soft direction. 

With the above justification, we can employ a some- 
what less involved path integral formulation us- 
ing the following boundary conditions: i5n^ • x^ — 0, 
dzSn^ ■ y^ = 0. The partition function of the fluctuat- 
ing director field subject to the above constraints can be 
written as 

Z = J VSn{r) 5{5n- ■ x-}S{dzSn- ■ y"} (4) 
X 6{6n+ ■ ±+}S{d,6n+ ■ A+} e-^/*^^^. 

The functional delta functions can be written as integral 
representations by introducing four Lagrange multiplier 
surface fields: 



Vr^ [h{r^) Sn-{r^) ■ x" + Ur^) djn-{r^_) ■ y- + Z3(r^) Sn+{v^) ■ x+ + k(r^) djn+(r^) • y+] | , 



where k = K/ksT. The integration over the director field is now Gaussian and can be easily performed. It yields 
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where is the angle between the corresponding soft axes 
of the two plates (see Fig. 1). The remaining integra- 
tion over the Lagrange multiplier fields can be performed 
to give InZ = — i X]q^ (det M(qj_)), which leads to a 
simple expression for the free energy per unit area: 
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Integration over u then gives the final result Eq. (|l|) 
above. The function f{9) = X^fcli cos'''(20)/fc'^ which 
describes the orientational dependence of the interaction 
has the structure of a zeta function that is commonly 
present in "Casimir" interactions. 

It is instructive to examine the limiting cases of plates 
in which the corresponding hard and soft axes are parallel 
or perpendicular to each other. For 9 — the boundary 
conditions corresponds to a hard-hard configuration for 
one component of Sn and a soft-soft one for the other. 
One can check that /(O) = ((3) ~ 1.20206 (C is the zeta 
function), thus we obtain exactly the same expression for 
the "Casimir" energy as in Ref. for "alike" boundary 
conditions. On the other hand /(■!) = — |C(3), which 
gives the same result as in Ref. for "unlike" boundary 
conditions {9 = 7r/2 corresponds to a soft-hard configu- 
ration for both components of Sn). 

Our calculation suggests two kind of experiments: (i) a 
direct measure of the torque exerted between two plates 
at a fixed separation, and (ii) a measure of force as a 
function of separation for plates at various angles. 

A possible experimental setup for the direct observa- 
tion of the fiuctuation-induced torque could be a torsion 
pendulum similar to the one discussed in Ref. |2^ ]. Our 
results imply that the torsion coefficient of the pendu- 
lum k (defined as the ratio between the torque t and the 
angular rotation A9) is corrected at = by an amount 



12 d2 



(9) 



due to the "Casimir" effect. Here R is the radius of 
the plates of area nR'^. Using the typical values knT — 



4.1 X 10"^^ dyncm, R — 1.5 cm, d — 10~^ cm, one ob- 
tains Afc ^ 10^^ dyncm. This accuracy may be reach- 
able using modern micro-manipulation techniques. 

A measure of the force-distance relation for various 
angles 9 could also be performed. An advantage of this 
procedure, as compared to measurement of the "simpler" 
effect corresponding to isotropic anchoring, is that rela- 
tive effects are more easily detectable. Indeed, while the 
weak signal of the "Casimir" force can be swamped by 
stronger effects (Van der Waals, etc.), the difference be- 
tween measurements performed at = and 9 = 90° 
should provide a differential evidence of the Casimir scal- 
ing. 
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